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In this note we address the problem of unconventional charmonium-like levels from the standpoint 
of level spacing theory. The level distribution of the newly discovered vector resonances is compared 
to that of standard charmonia analyzing their spectral rigidities. It is found that the unconventional 
charmonium-like states are significantly more compatible with the hypothesis of being levels from 
a Gaussian Orthogonal Ensamble of Random Matrices than the standard ones, which in turn seem 
more likely to be Poisson distributed. We discuss the consequences of this result and draw some 
hints for future investigations. 

Introduction. With the very recent observation of the charged resonances 2'(10610) and Z(10650) by 
the Belle collaboration [T], the family of unconventional quarkonium-hke states has further grown. Since the 
discovery of the X(3872), now observed also in LHC experiments, a long list of new narrow resonances has 
been found. There is a vast consensus that most of them are multiquark structures although a general picture 
is still missing. Some of these states occur extremely close and some other far from open-charm or beauty 
thresholds. For the close-to-thresholds ones, several authors agree that the appropriate interpretation is in 
terms of S"— wave 

£)(*)!)(*) or hadron molecules with a very small binding energy (compatible with 

zero), yet rather stable to be as narrow as the observation shows. Also the prompt production of AT (3872) in 
pp collisions at CDF has been observed making at least questionable the chances of a loosely bound molecule 
interpretation [5]. On the other hand it has been claimed that final state interactions mechanisms could be 
at the core of the surprising stability of such a molecular object [3]. Similarly the newly discovered states, 
the Z(10610) and Z(10650), have immediately been interpreted as hadron molecules ^ for their mass values 
happen to be exactly at the threshold values of BB* and B*B* mesons. 

Hadron molecules are meant to be extended tetraquark objects (several fermi in size) in which the strong 
interaction is conveyed by some long range pion exchange or rescattering mechanism. As opposed to this 
picture one could theorize the existence of compact tetraquark structures which are just new kind of hadrons 
with four quarks neutralizing the color within the typical range of strong interactions [S]. In principle, compact 
tetraquarks are not expected to be formed at the mass values of meson molecules; on the other hand some QqQq 
bound state could fluctuate into {Qq){Qq) or {QQ){qq) and the discrete levels of the unknown Hamiltonian 
binding QqQq should, as a result, be coupled to hadron molecule levels. 

In this letter we test the assumption that the known 1 resonances located at the mass values £ = 
{3943,4008,4263,4360,4634,4664} MeV - ah of them candidates to be exotic hadrons [1] - represent the dis- 
crete levels of some unknown compact tetraquark Hamiltonian along the same lines as the standard charmonia 
at 5 = {3096 (J/V'),3686 (i/'(2S')), 3772, 4039, 4153, 4421} MeV are the levels of the cc Hamiltonian with the 
Cornell potential. Resonances in £ are all produced in e~^e~ collisions with initial state radiation. Most of the 
levels of the exotic set £ happen to be away from open charm threshold and thus represent a good laboratory to 
explore the possibility that we are observing the spectrum of a complicated multiquark Hamiltonian; for earlier 
attempts of this kind see [7] . 

A very much studied conjecture in the field of quantum chaotic systems [B], states that spectra of quantum 
Hamiltonian systems whose classical analogs are described by (strongly) chaotic Hamiltonians show locally the 
same fluctuation properties as predicted by the so called Gaussian Orthogonal Ensemble (GOE) for large di- 
mensional Random Matrices. A portion of the quantum Hamiltonian spectrum is rescaled to spacing one and 
the levels so obtained turn out to be distributed as the eigenvalues of the GOE Random Matrices in the limit 
of large dimensions. In this limit indeed the local properties of Random Matrix eigenvalues are extracted, as 
the Wigner semicircle appears locally flat. As a consequence the probability distribution of the level spacings 
is expected to follow closely the Wigner law 



W{s) = 




The eigenvalues of GOE matrices following this distribution show the typical level repulsion features studied at 
length in the context of nuclear resonances. 

Naive formulations of tetraquark semiclassical 4-body Hamiltonians are possible, for example relying on 
one-gluon-exchange models. Most likely all of them express a chaotic classical dynamics. The Hamiltonian 



describing the cc system is, on the other hand, very close to an integrable one. Thus it is expected to have the 
level clustering features of the Poisson spacing distribution P{s) = exp(— s), or at least a discrepant behavior 
with respect to the GOE eigenvalues. 

Using the tool of the spectral rigidity also known as the A— Statistics developed initially by Dyson and Mehta 
we study the short sets £ and S in the attempt of confirming or disproving the picture according to which the 
£ levels should more markedly match the expected behavior for the GOE ones than standard charmonia, S, do. 
We surprisingly find that this is indeed the case although our explorative analysis has its natural limit in the 
very limited amount of data at hand - the method of A— Statistics has been systematically applied for example 
in the discussion of nuclear resonances level spacing where the data sets contain order of hundreds of levels. 

Yet we believe that this result is to be interpreted as an interesting suggestion which leads us to some 
speculative considerations we are still working on: i) exotic hadrons (for example those in the £ set) are 
just like the cc ones but with an additional light quark qq component; ii) they fall on the levels of some 
tetraquark Hamiltonian; Hi) once a discrete tetraquark level happens to be located within the level width of a 
molecular level - centered at some meson; threshold - because of the coupling between the two spectra induced 
by fluctuations like QqQq — {Qq){Qq) — > QiQq, the molecular level, otherwise very broad, gets metastable 
because of a Feshbach-like mechanism. 

^Statistics . The spectral rigidity (SR) is a measure of the deviation of a level set from uniform spacing: 
the more regular the set, the smaller the value of the spectral rigidity. Consider a set of N levels {Ei} rescaled 
to unit spacing, namely E^ — Ei = N — 1 = 2L, and centered with respect to the origin {Ei = —L and En = L). 
The sample cumulative function is 

c{x) = ^ e (x -£;,)- 5] e{E, - x) 

Ei>0 Ei<0 

where Q{x) = 1 if a; > and Q{x) = if a; < 0. The spectral rigidity, in its original form due to Dyson and 
Mehta, is defined as 

A3 ^ ^ min £^ (C(,s) -As- sf ds (2) 

following the notations introduced in [S] . 

The conjecture above means that a sequence of N experimental levels has to be compared with a sequence 
of N eigenvalues extracted from an ensemble of random matrices with large dimension D. Calculations in |10j 
show that in the large D limit the mean of the SR computed with the Poisson distribution is linear in the 
number of spacings, while that computed with the GOE ^13j eigenvalue distribution grows only logarithmically: 
it is therefore possible to discriminate between Poisson and GOE levels. This discrimination is more effective 
as the number of consecutive levels in the studied sequence grows. In practice, the number of levels available 
is often too low for A3 to provide a clear discrimination between Poisson and Gaussian Orthogonal Ensembles. 
It is therefore useful to consider a different notion of SR in order to reduce the variance. Following Bohigas et 
al. jTT] we set 

A3(a;, y) = - min / (C(s) -As- B? ds (3) 
y J.^ 

where A3 (a;,?/) is a generalization of the Dyson- Mehta estimator, recovered as y = 2L and x — —L. We thus 
define a new random variable A{w, y) built on averaging the spectral rigidity of smaller portions of the dataset 

pL+w-y 

where y takes continuous values between and 2L + 2w, which is the number of spacings of the sequence, while 
X ranges from —L — w to L + w — y. We have thus defined a family of statistical variables depending on y and 
w. A3 looks at the whole data set whereas A3(a::, y) checks a smaller number of levels and A{w, y) is his average 
on the data set. The parameter w is introduced in order to minimize possible finite size effects. The original 
definition of A3 is recovered in the limit w = , y 2L. 

Results. The data sets at hand are £ = {3943,4008,4263,4360,4634,4664} MeV, the candidate ex- 
otic levels, and the standard charmonia at 5 = {3096 (J/V'),3686 (V'(2S')), 3772, 4039, 4153, 4421} MeV. 
£ contains the masses in MeV of the l"" states G(3900) , ^(4008) , ^(4260) , r(4360) , X(4630) and 
y(4660). S contains the masses of the standard 1 charmonia, from the J/ip to the -0(4415). A third 
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set 8' = {3943,4008,4263,4360,4661} MeV is composed by the same resonances of the first set f , where the 
F(4630) and the F(4660) are taken to coincide with the Yb(4660) state, as proposed in [T7. In this study we 
will neglect the uncertainties on the masses. 

In order to choose the parameter w we are introducing, we study the behavior of A(i(j, y) on some test series. 
We choose w = 1, because smaller values are insensitive to variations of the spacings at the extrema of the 
series, whereas greater values are useless as they do not add further information. 

From now on we will use the notation A(l, y) = A(y). As we deal with 5 and 6 level sequences, we generate a 
large number of h.{y) samples from GOE and Poisson series. The GOE samples are obtained by diagonalizing 
30 random GOE matrices 4000 x 4000 in size, obtaining 24000 series of 5 levels each and 19980 series of 6 
levels each; the number of Poisson samples is similar (24000 and 20000). The integral in /^^{x^y) is evaluated 
analytically, whereas the A's are obtained by a midpoint rectangle approximation. We choose to evaluate A(?/) 
for integers and half-integers, with y £ [0.5, 7] (A(0) is identically zero). Note that each experimental data set 
has to be compared with samples of same cardinality. 

A comparison between the statistical properties of A3 as a function of the number of spacings and our A{y) is 
given in Figure IT] There we introduce the ensemble averages (•)p qqe where P stands for averaging against the 
Poisson Ensemble and GOE is for the Random Matrix ensemble. The ensemble average (A3)p ^^-^ is computed 
exactly in [10]. The (A(j/))p ^^^-^ is computed via the Monte Carlo sampling described above. 

The (A(?/))p have respectively the same linear and logarithmic behavior as (A3)p ^q^,, apart from the last 
points where finite-size effects dominate. A clear discrimination between GOE and Poisson sets is reached at y 
large. 
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FIG. 1: The solid lines represent {A3)p qq^, for series of y + 1 levels. The points show the Monte Carlo sampled 
(A(y))p,GOE- 

We now study the level properties of the experimental series £ and S. We observe that it is hard to discriminate 
between Poisson and non-Poisson sets because of the large variance of the Poisson A{y) random variable. On 
the other hand the GOE distributions have a smaller variance, so a more significant discrimination is possible 
between GOE and non-GOE sets by looking at high values of y. 

In Fig. [2] we show the experimental A(y) for the £ and S levels compared to the six level GOE averages. The 
results obtained for the sets £ and £' are compatible with the hypothesis of GOE distributed levels whereas 
this turns out not to be true for the S set. 

Given the random variable A(y) relative to the GOE ensemble we compute the related distribution function 
fy as depicted in Fig. [3[ We thus introduce 

ay = / fy{s)ds (5) 

J{s:fy{s)<fy} 

where fy is the value assumed by the distribution fy in correspondence of the SR A{y) associated to the 
experimental data set S. ay takes values in the [0, 1] interval; in correspondence of small ay values the null- 
hypothesis that 5 is a realization of the GOE ensemble has to be rejected. 
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FIG. 2: In the two figures the experimental and sampled data for 5-levels (left panel) and 6-levels series are shown. 
Solid lines show the experimental A(y) for the 5, E and E' series. For any j/, the black line indicates the point where 
the maximum of the distribution jy is attained; the colored areas indicate the regions that contain 50% or 90% of the 
A(y) samples. Note that when A(y) is outside the 90% area, the parameter a.y is smaller than 0.1. 



We consider the most significant five cases y — 5,5.5,6,6.5,7 and compute the corresponding cty by con- 
structing a binned distribution jy. By averaging over different binning choices we obtain the data reported in 
Table. 1 where for y >5 the fraction ay is smaller than 0.1. 



y 


5 


5.5 


6 


6.5 


7 


ay 


0.097 


0.084 


0.095 


0.096 


0.089 



TABLE I: ay is the fraction of sampled A(j/) whose probability is less than the value / introduced in the text. 



Computing the analogous quantity for the data sets E and £' we get values larger than 0.4 for y > 5. 

It is then reasonable to reject the hypothesis that S levels are extracted from the Gaussian Orthogonal 
Ensemble. Note that, because of the large variance of the Poisson samples, it is practically impossible to reject 
the hypothesis that a series of few levels (< 20) belongs to the Poisson Ensemble, but, as observed above, it is 
easier to state that they are not of the GOE type. 

As an example in Fig. [s] the sampled distributions for A(7) are shown and the experimental values of the S 
and E sequences. 
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Poisson samples 
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A(7) Exotica 



A(7) 



FIG. 3: Sampled distributions for A(7). The histograms are obtained with 19980 GOE samples and 20000 Poisson 
samples and rescaled to unit area. The vertical lines indicate the A(7) value for the standard and the exotica series. 

The bin width is equal 0.008. 



Conclusions . Comparing the spectral rigidity of standard 5'- wave charmonia with that of the unconventional 
vector resonances recently discovered, we observe that the latter states are more significantly compatible with 
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the hypothesis of being the levels of some multiquark Hamiltonian (whose classical analog exhibits chaotic 
dynamics and therefore having quantum levels distributed with the Wigner law) than the former, which, on 
the other hand, could be thought as the levels of some classically integrable one - as the simplest version of 
the Cornell potential Hamiltonian is. The limit of this analysis is in the small amount of experimental data 
available both in the standard and unconventional sectors. We have introduced a slight modification of the 
spectral rigidity estimators used in the literature to improve as much as possible the quality of our analysis with 
a small number of levels. 

Molecular Hamiltonians, besides the fact that describe two-body systems, could as well be regulated by 
complicated potentials with Wigner distributed quantum level spacings. Yet there are states in the £ sequence 
of unconventional 1 resonances which do not match molecular thresholds whereas the most accredited hadron 
molecule model has S'-wave molecules almost exactly at threshold. There are no clear hints on the form of these 
potentials neither and the main problem of the spectroscopy of the new X, Y, Z resonances remains that of 
finding a unified description that accounts for both on and off-threshold particles. 
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